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Based on the exact solution of the eigenvalue problem for the Z7q [s/ (2 1 1)] vertex model 
built from alternating 3-dimensional fundamental and dual representations by means of the 
algebraic Bethe ansatz we investigate the ground state and low energy excitations of the 
corresponding mixed superspin chain for deformation parameter q = exp(— ij/2). The model 
has a line of critical points with central charge c = and continua of conformal dimensions 
grouped into sectors with 7-dependent lower edges for < 7 < ir/2. The finite size scaling 
behaviour is consistent with a low energy effective theory consisting of one compact and one 
non-compact bosonic degree of freedom. In the 'ferromagnetic' regime n < 7 < 2tt the critical 
theory has c = — 1 with exponents varying continuously with the deformation parameter. 
Spin and charge degrees of freedom are separated in the finite size spectrum which coincides 
with that of the U q [osp(2\2)] spin chain. In the intermediate regime tt/2 < 7 < ir the finite 
size scaling of the ground state energy depends on the deformation parameter. 

I. INTRODUCTION 

Studies of exactly solvable two-dimensional vertex models or the equivalent (1 + l)-dimensional 
quantum spin chains can provide important insights into the nature of excitations in strongly 
correlated systems and their critical behaviour. Over the years this approach has provided much 
to the present understanding of such models based on ordinary Lie algebras whose massless regime 
is believed to be described by Wess-Zumino-Witten models on the corresponding group. On the 
other hand, many of the physical properties of vertex models based on Lie superalgebr&s and 
their quantum deformations are still not understood in detail. At the same time and in spite of 
significant progress in recent years only little is known about the likely candidates for the low energy 
effective description of these vertex models, i.e. (1 + l)-dimensional conformal field theories with 
non-compact target spaces. Advances in this direction are highly desirable as they would likely 
lead to progress for some problems related to the duality between gauge and string theories (see 
Ref. Q and References therein) but also in statistical mechanics, e.g. for the description of disorder 
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driven phase transitions within the superspin approach to non- interacting electron systems, see 
e.g. fl 3- 

One possible approach to this problem is based on the observation that a non-compact con- 
tinuum limit can arise from lattice models with a finite number of states per site [j-6]. If such 
models are integrable the powerful techniques of the quantum inverse scattering method allow for 
a detailed analysis of their spectrum and ultimately provide important insights into their contin- 
uum limit. Concerning the possible applications mentioned above lattice models with alternation 
between conjugate representations of the superalgebra have been found to be particular important. 
For the integrable sZ(2|l) superspin chain mixing the fundamental representation 3 and its dual 3 
this approach has led to the identification of the continuum limit with the SU (2\1) WZW model 
at level k = 1 [J]. Recently, this has been generalized to find the scattering theory arising in the 
continuum limit of the antiferromagnetic gl(n + N\N) spin chains with n, iV > [7]. 

Another direction in which these results may be generalized is by deformation of the underlying 
symmetry: in the case of ordinary Lie algebras this has led to models which exhibit critical lines 
with anomalous exponents depending continuously on the deformation parameter. Whether such 
a behaviour occurs in superspin chains when a deformation parameter is introduced and even 
whether the critical behaviour of the mixed superspin chain observed in the undeformed case is 
robust against the deformation is the question which we want to address in this paper. 

Our paper is organized as follows: below we recall the definition of the mixed C/ g [s/(2|l)] vertex 
model [8] from which we obtain the integrable superspin Hamiltonian which is then solved by 
means of the algebraic Bethe ansatz. Since the analysis of the Bethe equations makes use of the 
known properties of the Fateev Zamolodchikov model [9] with twisted boundary conditions we also 
summarize what is known about the different critical phases of the latter. In Section HVl we present 
our results on the low energy properties of the mixed superspin chain in the 'antiferromagnetic' 
regime which can be interpreted as a regularization of a continuum theory consisting of a compact 
and a non-compact boson, similarly to the sZ(2|l) mixed superspin chain discussed in Ref. |4. 
The critical behaviour in the 'ferromagnetic' regime discussed in Section [V] turns out to be very 
different: the corresponding low energy theory exhibits exact spin-charge separation for all values 
of the deformation parameter q. Both the spinon and the holon modes are described by U(l) 
Gaussian theories. In the isotropic limit restoring s/(2|l) invariance the spin part of the spectrum 
acquires a quadratic dispersion above an (completely polarized) reference state while the charge 
part of the s pec trum remains conformal in the same universality class as the isotropic osp(2\2) spin 
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II. MIXED VERTEX MODEL 

The weights of the mixed vertex model are based on the i?-matrices associated to the three 
dimensional vector representation of f7 g [sZ(2|l)] and its dual, labelled 3 and 3 in the following. 
These i?-matrices act on the tensor products 3® 3, 3 <8> 3, 3 <8> 3 and 30 3 and we shall denote them 
by R( 33 \X), R( 33 \A), BS 3 ' 3 \\) and i?( 3 ' 3 )(A) respectively. For a general discussion of i?-matrices 



alternating between the vector representation of U q [sl(n\m)] and its dual see for instance [12M15I]. 
In the specific case of the quantum superalgebra ?7 9 [sZ(2|l)] the above set of i?-matrices have been 
explicitly discussed in Q] for a particul ar g rading (models of this type without grading have been 



introduced before by Perk and Schultz 



161]). In what follows we shall present their expressions for 



arbitrary ordering of the Grassmann parities, 
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<?(A) = £ «,(A)eg eg + 6(A) £ eg eg + c(A) { 
j=i j,k=i 



E(-» 



> W '«S ® eg 



+ exp(-2A) £ (-l)* Pfc e$ ® e S |> . (2-4) 

where eg G End(C 3 ) are the standard 3x3 Weyl matrices. The symbol pj denote the Grassmann 
parities distinguishing the bosonic pj = and fermionic pj = 1 degrees of freedom. 

The dependence of the Boltzmann weights Oj(A), 6(A) and c(A) on the spectral parameter are, 

sinh [A - i(2pj - 1)7)] sinh [A] sinh [i 7 ] 

^ (A) = sinh[A + * 7 ] ' b(A) = sinh[A + n ] ' ^ = ^ [A + i 7 ] ' M 

while functions /(j, A;) depend only on the anisotropy 7 as follows, 

/(l,2)=exp[2*7(l-p 3 )], /(l,3)=exp[2z7p 2 ], /(2, 3) = exp[-2i 7 pi]. (2.6) 

The i?-matrices defined above fulfill the Yang-Baxter equation on any tensor product built up 
from the 3 and 3 representation, namely 

R[^\\)R{f^\x + ^a^V) = 4^ 3 V)43 1,W3) (A + /^'^(A), (2.7) 

where the representations CJj £ {3, 3} for j = 1, 2, 3. 

One consequence of these set of Yang-Baxter relations is that there exists two different types of 
Lax operators obeying the Yang-Baxter algebra with the same i?-matrix. For example, this means 
that an integrable vertex model combining the i?( 3 ' 3 ) (A) and i?( 3,3 ) (A) can be constructed within the 
framework of the quantum inverse scattering method. As usual the respective row-to-row transfer 
matrix is written as the supertrace [l^] over the auxiliary space A ~ C 3 of the following ordered 



product of operators: 



T( 3 )(A, = Stv A ^' 3 2(A)^ ) -i(A - »7 + (A) • • • ^2i 3) (A - H + 



(2.i 



acting on the Hilbert space (3(8) 3)® ~ C 2L . The alternation on the spectral parameter can be 
introduced since the i?-matrices are additive on A. Note that this choice of inhomogeneity does not 
spoil the basic properties such as the symmetry and locality of the interactions of the corresponding 
alternating superspin chain with the Hamiltonian obtained by taking the logarithmic derivative of 
T( 3 )(A,£) at A = 0. 
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By the same token a solvable integrable vertex model alternating the /^-matrices i?( 3 ' 3 ) (A) and 
i?( 3 ' 3 )(A) can be constructed. The expression of the transfer matrix mixing such operators is 



(A, = Str A K^l (A + i 7 - O^S-i (A)^S_ 2 ( A + i 7 - • • • ^ (A) 



(3,3) 



,(3,3) 



,(3,3), 



(2.9) 



acting on the same Hilbert space as (|2.8p above. Again, an alternating superspin chain can be 
constructed by expansion of the transfer matrix around A = 0. 

It turns out that - in addition to commuting among themselves - the transfer matrices (A, £) 
and T( 3 )(A,0 constitute a family of commuting operators when the inhomogeneities £ and £ are 
the same, i.e. 



0, VA,/i. 



(2.10) 



The property (|2.10|) relies on the fact that we have chosen an identical ordering of representations 
3 and 3 in the definition of the Hilbert spaces for the transfer matrices (|2.8|) and (|2.9|) and follows 
from the Yang-Baxter equation (|2.7p once we choose the representations uj\ = 3, uji = 3 and 
= 3,3. In this situation we are able to construct an integrable vertex model that alternates 
representations 3 and 3 both on horizontal and vertical spaces of states of a square lattice of size 
2L x 2L. The respective 'double row' transfer matrix of such model is obtained by taking the 
following product, 



T {mix) q = T (3) ( A) ^ r (3) ( A) 



(2.11) 



By construction T^ mix \\ = 0,£) is proportional to the translation operator by two lattice sites. 
Therefore, we can define an integrable superspin Hamiltonian, 



H (mix) = • lnT (™z)( A) £) 
OA 



(2.12) 



A=0 



The expression for Jj( mix ) ' m terms of the /^-matrices can be obtained by computing the individual 
Hamiltonians associated with the transfer matrices T"( 3 )(A) and T^ 3 )(A). The technicalities entering 
this computation are cumbersome but the final result is somehow simple, 

2L 



3=2 

mod 2 
2L 



rcg&K - i 7 ) + p jd+2 n^ 2 mm^ ~ <7) 



(3,3) 



,(3,3) 



+ I 



nf^xiil -0+ Pjj+iKff^myf/Sl - (2-13) 



,(3,3) 



(3,3) 



3=1 

mod 2 
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where periodic boundary conditions 2L + 1 = 1 and 2L + 2 = 2 are assumed. The operator P a b 

3 



is the graded permutation P c 



ab 



j,k=l 



("I) 



PjPk 



(a) 
I k 



e^ 1 - and 1Z a i,(\) denotes the derivative of the 



i?-matrix lZ a b{\) with respect to the spectral parameter A. 

The diagonalization of the above transfer matrix can be carried out by applying the nested 
algebraic Bethe ansatz approach [lsl [lij]. For this particular mixed vertex model the essential 
tools to obtain the eigenvalues of T^ mix \\, £) can for instance be found in [20J. We shall not repeat 
here these technical details and concentrate our attention only to the main results. As usual 
the expressions for the eigenvalues obtained in this approach will depend on the choice of grading 



13, 



21 



24 ] . For later convenience we use [pi,P2,P3] = [0, 1, 0] in the following. Let A^™^ 2 (A) denote 



the eigenvalues of T^ mix \\, £) in the sector of the Hilbert space selected by fixing the two conserved 
quantum numbers related to the U(l) subalgebras of U q [sl(2\l), i.e. charge b = (iVi — A r 2)/2 and 
z-component of the spin S3 = L — (N\ + A r 2)/2. As a consequence of ()2.10j) the eigenvalues can be 
factorized 



aSW=a?;, 2 (a)a^(a). 



,(3) 



(3) 



(2.14) 



Here A N N (A) and A N N (A) are the corresponding eigenvalues associated to the transfer matrices 
T( 3 )(A,£) and T®(A,£) respectively. 



It turns out that the expressions for the eigenvalues A Ni N2 (X) and A^ N ^) (A) are given by, 



,(3) 



A 



(3) 

N lt N 2 



(A) 



and 
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+ 



A - i 7 /2) 
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sinh(A) ] L -^-sinh(A-Af ) +i 7 ) 



sinh(A + i'y) 
- »7/2) 



g sinh(Af 



sinh(A- Xf ] ) 
- A - i'y) 



1= \ sinh(A - Xf ] + i 7 /2) fJi sinh(Af } - A) 



(2.15) 
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(3) 
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(A) 
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sinh(A) 
sinh(A + i'y) 



L Q sinh(A - 
j =1 sinh(A 



A« 



£ + 3*7/2) 
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sinh(A — £ + i'y) 
sinh(A — £) 



L n 2 

n 



sinh(Af } — A + £) 



sinh(A 



(2) 



A + £ - * 7 ) 



sinh(A — £ + ij) sinh(A) 
sinh(A + i'y) sinh(A — £) 



L Nl sinh(AS 1) 



sinh(A 
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where the rapidities Xj and satisfy the following set of nested Bethe equations, 



sinh(Aj 1) + *7/2) 
sinh^ - ij/2) 

sinh(Af } + £) 

(2) 

sinh(A!- + £ — 27) 



* sinh(Aj 1) 
11 sinh(A« 



*£° + *7/2) 



fe=l 



A? - ry/2) 



v(2) 
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^ sinh(A} 

M 8inh(Af - Ap - z 7 /2 
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(2.17) 



Note that for the particular choice £ = « 7 /2 these Bethe ansatz equations become symmetrical 
on the variables Xj and A*- . In the following we will concentrate our studies on this case. The 
eigenvalues of the Hamiltonian (|2.12p corresponding to a solution of (|2,17p are 



h N 1 ,N 2 \^) 







/ ATi 



4Lcot^ + 2(]T 



A=0 
sin 7 



sin 7 



N 2 

< <)S ") ~ ( os '' "^Ag 5^ COS7 — cosh2Ag 



(2.18) 



Considering the above solution one sees that the spectrum at the points 7 and 2tt — 7 are related 
to each other by only a sign, 



E (mix) , x 



-E 



(mix) 
N U N 2 



(27T - 7 ) . 



(2.19) 



III. THE TWISTED XXZ SPIN-1 MODEL 



The classical vertex model associated to the integrable Heisenberg XXZ spin-1 chain turns out 
to be the three-state factorized i?-matrix found by Zamolodchikov and Fateev [9J]. Considering our 
previous notation this operator can be expressed as 



n ab {x) 



=eg« 
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where the corresponding Boltzmann weights 6(A), c(A), d(X), /(A) and /i(A) are given by, 
sinh(A) sinh(^7) - sinh(i7) sinh(A) 



6(A) = 
/(A) 



elA) 



d(X) 



sinh(A + i7)' sinh(A + i 7 ) ' sinh(A + iZ) sinh(A + £7) ' 

sinh(A - q) sinh(A) _ 2cosh(f ) [sinh(i^)] 2 



sinh(A + iZ) sinh(A + 17) ' 



MA) 



sinh(A + «?) sinh(A + 17) 



(3.1) 



(3.2) 
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With (|3.ip the respective transfer-matrix T(X) with toroidal boundary conditions can formally 
be constructed as follows 



T(A) = Tr A [GaRalWR-al-iM • • -^i(A)] 



(3.3) 



where Ga denotes a 3 x 3 matrix representing the twisted boundary condition. 

The diagonal twisted boundary condition compatible with integrability is obtained by choosing 
the matrix G as, 

I 1 ^ 




e 



2i<p 



(3.4) 



where the angle (p is assumed to be in the interval < ip < ir. 

The transfer matrix (|3.3p with (|3.4[) can be diagonalized with very little difference from the 
periodic case since the presence of the diagonal boundary matrix Ga preserves the U(l) bulk 
symmetry. The respective eigenvalues can be determined either by using the mechanism of fusion 
2514271] or by applying the algebraic Bethe ansatz construction developed in |28j . As a consequence 
of the U(l) invariance of the transfer matrix the Hilbert space can be separated in disjoint sectors 
corresponding to total magnetization S3. Starting from the state with maximal S3 = L one obtains 
the following expression of the corresponding eigenvalues A^(X,ip) in the sector S3 = L — N , 
N = (),••• ,L 



N 



J-Jj sinh(Aj — A — if/ 2) 



sinh(A — if/2) sinh(A) 
sinh(A + if) sinh(A + if/2) 



L N 

n 

5=1 



+ e 



•up 



sinh(A) 
sinh(A + if) 



L N 

n 

j=i 



sinh(A — Xj + if) sinh(A — Xj — if/2) 
sinh(A — Xj + if/2) sinh(A — Xj) 



sinh(A — Xj + if) 
sinh(A — Xj) 



(3.5) 



where the rapidities Xj satisfy the following Bethe ansatz equations, 



sinh(Aj +if/2) 
sinh(Aj — if/2) 



n 

k=l 



sinh(Aj — Xk + if/2) 
sinh(Aj — Xk — if/2) ' 



;" = !,••• ,N. 



(3.6) 



At this point we have gathered the basic ingredients allowing to establish a mapping among 
part of the spectrum of the mixed transfer matrix ()2.1ip and the eigenvalues of the XXZ spin-1 
model with special toroidal boundary condition for the special choice of £ = if/2 in the transfer 
matrix 

T (mix)^ of the mixed 

spin chain: as mentioned above the Bethe ansatz equations (|2,17p 
become symmetric between the two levels for this choice of £ and their structure resembles that of 
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the Bethe ansatz equations (|3.6|) of the XXZ — 1 chain with twist ip = ir. In particular we find 
that in the sector N\ = N2 = N of the Hilbert space of the mixed chain (this is where the total 
t/g[sZ(2|l)] charge b of the Bethe state is zero) there exists a subset of eigenstates parametrized by 
Bethe roots which can be identified with eigenstates of the spin-1 chain in the sector S3 = L — N 
by setting A^ = A,- = A,-. A similar correspondence has been observed between a subset of 
eigenvalues of the alternating si(2|l) superspin chain obtained in the limit 7 — > from the mixed 



chain considered here and the ST/(2)-invariant spin-1 Takhtajan-Babujian chain (4]. For such 
states, a direct inspection of the expressions for the eigenvalues (|2.15p . (j2. 16|) . ()3.5p leads us to the 
following relation, 

2L 



n.n y A > 



sinh(A + i7/2) 
sinh(A — ij/2) 



[A N (X, i p = 7r)]\ (3.7) 



As a consequence, we obtain from (|2.18p a relation between the energy eigenvalues of the XXZ 
spin-1 chain with L sites and those of the J7 ? [si(2|l)] superspin chain with 2L sites of alternating 
representations 3 and 3: 

E^f = 4L cot( 7 /2) + 2E$ XZ) (99 = vr) (3.8) 



where 



d 



E N XXZ \cp)=i ^ln A n (A, y ) 



y — 

^ — J COS 'y — rrwan 9 A j 
k=l 



A=0 7. 77, cos 7 ~ cosn 2 ^fc 



is the energy eigenvalue of the XXZ spin-1 chain corresponding to a solution of the Bethe equations 
(|3.6p . As in Eq. (|2.19p for the mixed superspin chain the spectrum of the XXZ model is inverted 
at anisotropy 7 = tt, i.e. spec( 7 ) — spec(27r — 7). 

In the thermodynamic limit L — > 00 the solutions of the Bethe equations (|3.6p are grouped into 



Am) 

a parity v m = dbl: 



'strings' consisting of m complex rapidities \j characterized by a common real center and 



A (m) = A(m) +i l {m+1 _ 2j)+i _ {l _ Vmh j = 1, . . . ,m . (3.10) 



26, 



29]. Here we 



The allowed values of (m, v m ) depend on the anisotropy 7 in an involved way 
shall not go into the details of the string classification: for the range of anisotropies < 7 < ir 
that we are considering in this paper it is found that the most relevant root configurations solving 
can be organized into strings (1,+), (1, — ) and (2,+). 
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A. The disordered antiferromagnet regime of the spin-1 chain 



Most importantly, the ground state of the system without twist, 



lattice is a condensate of L/2 (2, +)-strings in this regime 25l. l26l. l30|| . In the thermodynamic limit 



if = 0, and on an even length 



one can compute the energy per site giving €00(7) = —2 cot (7/2). The finite size spectrum of the 



31 



32]: in the entire interval < 7 < ir 



XXZ spin-1 model without twist has been investigated in 
the spectrum has gapless excitations with Fermi velocity vf = 2tt/^. The central charge of the 
conformal field theory describing the low energy sector is c = 3/2, hence the ground state energy 
for even L scales as 

E^\<p = 0) - LeUl) = c + o(i) = + oil) . (3.11) 

The operators of this CFT are given by products of Ising operators and U(l) Kac-Moody fields. 
The scaling dimensions of these composite fields in the presence of a twist (p are [33j 



Xj(0,0) €{0,1}, X / (0,1)=X J (1,0) = ^, X/(l,l) = i. 

Depending on the parity of L the possible subset of the KM representations is determined by the 
selection rules 

n = r + L mod 2, m = j + L mod 2 (3.13) 

for given parity j and toroidal b.c. (type r) of the Ising sector. The smallest exponents obtained 
from (|3.12p for even L are 

4 ; r /7r) (7)-^ 2 — 



'(0.0) w; "W 16X C ' 



<„r / ' , (7)=l+^+^^ — 



and for L odd 



^d.o) w; ~8 T CT W^ 16 X, 



X (i-0) 8 + ^1" -J YoX' 



As the twist is varied the energy of the tp = ground state ((r, j) = (0,0), (m,n) = (0,0) for 
even L) increases until there occurs a crossing with the level evolving from (n,m) = (0,-1) at 
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0.4 : 



> 0.2, 



-■ CFT 

♦ (n,m)=(0,0) 

• (0,-1) 
■ d,0) 

A (1,-1) L odd 
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FIG. 1. Evolution of the ground state and lowest excitation energies of the XXZ spin-1 model with 7 = 
as a function of twist if. Filled (open) symbols in black are for L = 6 (30), respectively; red symbols for 
L = 7 (31). Dash-dotted lines indicate the CFT predictions (|3.12p . 



(p = (3tt — 7)/4 (see Fig. [I]). The conformal dimension of this state at twist 92 = tt is Xq'® = 1/8 for 
even and X®'q =1/8 for odd length lattices independent of the deformation parameter 7. Together 
with the finite size scaling of the ground state energy (|3.1ip this gives a state whose energy is 
-^00(7) without any finite size corrections] Note that the corresponding wave function does vary 
with 7. According to (|3.9p this can only be realized with a highly degenerate configuration of 
Bethe roots, namely = for all k = 1,...,L. The identification (|3.8p implies the existence 
of a zero energy eigenstate of the mixed Z7g[si(2|l)] superspin model in the singlet sector of the 
latter. For < 7 < tt/2 this is the ground state of the 7r-twisted XXZ spin-1 chain and the mixed 
superspin chain, implying that the effective central charges of this models are c = 0. For 7 > tt/2, 
the two-fold degenerate level with scaling dimension X^' realized in the XXZ spin-1 chain of 

odd length has a lower energy at twist ip = tt . The consequences of this crossing for the superspin 
chain will be discussed below. 



B. The disordered ferromagnetic regime of the spin-1 chain 

As a consequence of the inversion of the spectrum under 7 -H- 7 = 2ir — 7 we can discuss the 
properties of the spin-1 chain in this regime in the same interval < 7 < tt while changing the 
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sign of the energies (|3.9|) . This leaves the string classification (|3.10p unchanged. Without twist 
in the boundary conditions the configuration of Bethe roots corresponding to the ground state in 
the disordered ferromagnetic regime is given by a filled sea of L (1, — )-strings. Above this state 
there are gapless low energy excitations with Fermi velocity vp = 2tt/(2tt — 7). The corresponding 
conformal field theory has been identified as a U(l) Gaussian model with central charge c = 1 and 
scaling dimensions of the primary operators given by [34J ] 

X n ,m{l) = n 2 x p + - — , x p = — (3.16) 

where n and m take integer values which determine the magnetization S3 = n and vorticity of 
the corresponding state. Note that the compactification radius of the boson R = ^Jxp~ vanishes as 
7 — > indicating the transition into the non-conformal isotropic ferromagnetic state. For toroidal 
boundary conditions with twist ip one has to replace m — > m + <p/2ir in Eq. (|3.16p . The adiabatic 
evolution of the Bethe roots under the twist is rather involved, a detailed study of the corresponding 
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regime in the spin-1/2 XXZ chain can be found in Ref. 

Here our focus is on the antiperiodic twisted chain: choosing 93 = it implies that the finite 

size gaps are given by Eq. (|3.16p with integer n but half-odd integer m = ±1/2, ±3/2, .... As an 

immediate consequence the lowest state of the conformal part of the spectrum in the finite system 
(xx y\ 

has an energy (ioo is the bulk ground state energy density of the chain in this regime) 

pm t~Axxz)^ nv F 2ttvf~ _ ttv f 2ttv f it 

Eo(L)-Le^ >~-_ + _X 0|J --— + — - (3.17) 

which grows as 7 — > 0, eventually leaving the range of applicability of the finite size analysis based 
on Eq. (gJBD . 

For additional insights into the properties of the spin-1 chain we have to rely on the analysis 
of the Bethe equations (|3.6p : the configuration of Bethe roots for this state evolves as 7 decreases 
from 7r to 0: for 7 G (tt/(& + l),n/k) we find that it consists of L — k (1, — )-strings and a single 
(k, +)-string according to the classification (|3.10p . see Fig.[2j In the finite system the configuration 
reduces to a single (L, +)-string for 7 < tt/L. This is a bound state of magnon-excitations over the 
ferromagnetic pseudo vacuum state with S3 = L. Based on this observation we propose that for a 
system of finite size L there is a level crossing at 7 = tt/L leading to a fully polarized ground state 
at smaller values of the anisotropy. This proposal has been confirmed by numerical diagonalization 
of the Hamiltonian. 
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Y=2jc/3 



y=2n/5 



y=2n/7 



-Till 



2 

Re(X) 



FIG. 2. Configuration of Bethe roots for the lowest S3 = state in the ferromagnetic disordered regime of 
the XXZ spin-1 chain with L = 8 sites for various values of 7. 

IV. ANTIFERROMAGNETIC REGIME OF THE MIXED CHAIN 

To discuss the properties of the mixed superspin chain we also need to distinguish two cases: 
< 7 < 7r: in analogy to the XXZ spin-1 chain we will call this regime 'antiferromagnetic', 
7r < 7 < 2tt: this 'ferromagnetic' regime will be discussed in the next section. 



The solutions to the Bethe equations for the mixed superspin chain (|2.17|) can be classified into 
strings - similarly as for the XXZ spin-1 chain in (|3.10p : in the large L limit Bethe roots with 
Im(Xj) 7^ or it/2 have to be combined such that their differences coincide with poles or zeroes 
of the bare scattering phase shifts on the right hand sides of Eqs. (|2.17|) . As in the rational case of 
the s/(2|l) mixed superspin chain [4] this considerations lead to 



'1) 1-strings 

just as for the XXZ spin-1 chain there are two types of unpaired roots allowed, namely real 
ones and roots on the line Im(X) = n/2, 



and to composites combining roots from both levels: 

(2) wide and strange strings 

these configurations consist of m roots from both levels of the Bethe ansatz with the same 
center X m G R U (R + in/2). For m odd they have been called wide strings in Ref. J, e.g. 



A 



A 



(i) 

m,k 

(2) 
m,j 



7 771 + 1 

X m + ij (m + 3 - 4k) , k = 1, . . . , — - — , 

7 m — 1 

X m + i-(m + I - 4j) , j = 1, . . . , — - — 



(4.1) 
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and a second type obtained by interchanging first and second level roots, {A^ fc } f-> {A m j}. 
For m even, the so-called strange strings have the unusual property of not being invariant 
under complex conjugation, e.g. 

^\ = X m + d(m + 3-4k), k = l,...,^, 

4 2 (4.2) 

= A m + iJ(m + l-4j), j = l,...,^. 

Again, there is a second type of such configurations for given m obtained by {A^ fe } f-> {A^ •}. 
(3) narrow strings 

Finally, there are composites which contain the same number m/2 of roots on either level. 
They may be seen as degenerations of two wide or strange strings with the same center A m . 

The existence of strings of a given length as well as their parity (i.e. whether they are centered 

around the real axis or the line Im(\ m ) = vr/2) depends on the deformation parameter 7, again as 

in the XXZ spin-1 chain. For the states with lowest (and highest) energies, we find that 1-strings of 

either parity and strange 2-strings centered around the real axis plus their possible degenerations 

into narrow ones are sufficient to capture the spectrum (in the 'ferromagnetic' high energy regime 

this is true at least for 7 > ir/3 as discussed in Section IVl below) . 

In the antiferromagnetic regime this classification is particularly useful in the sector where 

the total number of roots on the two levels is the same, i.e. for N\ = N2 in f|2. lTj) : here the total 

C/ g [sZ(2|l)] charge of the Bethe state is zero. Many low lying excited states in this sector correspond 

to root configurations consisting of strange 2-strings (up to finite size corrections), i.e. sets {A^} 
(2) 

and {A^ } that are mapped onto each other by complex conjugation. This sector also contains the 
lowest state in the singlet sector of the model, where N± = N2 = L. As has been argued above, 
this state is also present in spectrum of the XXZ spin-1 chain related by (|3.8p and has energy 
jjj(mix) = This is the ground state of the superspin chain for 7 < it/ 2. The corresponding root 
configuration solving the Bethe equations (|2.17p is a£ = for all k = 1, . . . , L and a = 1, 2. The 
same observation has been made in the rational model obtained as 7 — > [4]. 

A. Small systems 

For L up to 4 we have computed the complete spectrum of the mixed chain by exact numerical 
diagonalization of the Hamiltonian. As a consequence of the deformation some of the degenerations 
present in the s/(2|l)-symmetric superspin chain are lifted. For example, an sZ(2| 1) octet [b, s] = 
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TABLE I. Low energy states of the L = 3 superspin chain in the antiferromagnetic regime for 7 = 27r/7 
and the identified Bcthc configurations. Additional root configurations to these energies can be obtained by 
using the symmetries of the Bethe equations. 



{Ni,N 2 ) 


Energy E 


degeneracy 


Bethe roots 


(3,3) 





1 


A (i) 


= {0,0,0} = A< 2 ) (XXZ) 




(2,2) 


1.2968 


2 


A (i) 


= {±^0.2386} = A( 2 ) (XXZ) 




(3,3) 


1.6523 


2 


AW 


= {-0.0110 ± i0.2348, 00} = -A< 2 > 




(3,2) 


2.3027 


4 


A(D 


= {±i0.2492,i7r/2}, ^ (2) = {±i0.2255} 




(2,2) 


4.9748 


4 


AM 


= {±0.1161 + i0.2704} = (A< 2 ))* 




(3,3) 


5.1859 ±i0.6690 


2*2 


A(D 


= {-0.1429 - i0.2668, 0.0932 - i0.2675, 00} = 


_ A ( 2 ) 


(3,2) 


6.1084 ± i0.3325 


2*8 


A(D 


= {-0.0889 + i0.2822, 0.1199 + i0.3044, 0.2581 


- i0.8296} 








A (2) 


= {-0.0790 - i0.2295, 0.1415 - iO.2039} 





[0, 1] splits into two charge doublets with S3 = ±1 and S3 = respectively and a quartet with 
charge ±i and S3 



2" 



[0, 1] ^{\b = 0, s 3 = 1), \b = 0, s 3 = -1)} U {|6 = 0, s 3 = 0), \b = 0, s 3 = 0)} 



U{|6 



1 



• S3 



2 ) ' |6= 2' S3 



1 



),\b 



1 



1 



1 



• S3 



)}• 



(4.3) 



2 2' 1 2 2' 1 2 2' 1 2 2' 

Interestingly, we observe cases where pairs of the S3 = doublets arising from degenerate octets in 
the isotropic case split further into pairs with complex conjugate eigenvalues. 

For L = 3 we have identified the low energy states in terms of their corresponding configuration 
of Bethe roots, see Table [J for the spectrum at anisotropy 7 = 27r/7 (the degeneracies found in 
the numerical solution can be reproduced by applying symmetry operations on the set of Bethe 
roots, i.e. A^ 1 ) o A^ 2 ) or A( a ) — A^ a ) for both a = 1,2, and by using the global symmetries of 
the mixed chain, i.e. reversal of all spins). Although they are strongly deformed in some cases, the 
string content of these configurations according to the classification given above can be identified. 

Based on the numerical and analytical data the following general picture for the lowest excita- 
tions emerges: 

As shown in Ref. the low-lying multiplets of the s/(2|l) symmetric chain (7 — > 0) are the singlet 
ground state with E = in the normalization used here followed by a single sZ(2|l) octet as the 
lowest excitation. Above these there are two more degenerate octets and a pair of degenerate 
indecomposables, each containing 8 states. Upon deformation the degeneracies of these excitations 
lifted as described above, see Fig. [3l Note that the S3 = 1 doublet arising from the lowest octet is 
part of the XXZ spin-1 subsector of the spectrum which according to Eq. (|3,15[) has a finite size 
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1 2 



y 

FIG. 3. Low energy part of the spectrum of the L = 3 chain for anisotropy < 7 < 7r: black lines and 
filled symbols denote the real eigenvalues, red lines and open symbols the real part of complex eigenvalues 
of the Hamiltonian. Full and dashed lines denote eigenvalues obtained by solution of the Bethe equations, 
see Table HI whereas the symbols are eigenvalues obtained by numerical diagonalization of the Hamiltonian 
for which the corresponding configuration of Bethe roots has not been identified. Dotted lines connecting 
numerical data are guides to the eye only. 

energy gap AE^ XXZ ^ = (2-kvf/L) \X c — |] in the large L limit. For ir/2 < 7 < tt this doublet is 
the ground state of the mixed chain, as expected from the analysis of the XXZ-1 chain with anti- 
periodic boundary conditions for odd L. As 7 — > tt all states evolving from this octet degenerate 
at E = 0. 

The 16 states of the degenerate octets split into a quartet and an octet with real energies and 
two doublets with complex conjugate eigenvalues of the Hamiltonian. The latter eigenvalues also 
approach as 7 — > tt. Similarly the states of the degenerate indecomposables split into two quartets 
and an octet, all with real energies. Increasing 7 beyond ~ 7r/4 the real octets from these two groups 
degenerate and turn into two octets with complex conjugate energies. Our numerical data indicate 
that this conversion of pairs of real eigenvalues into pairs of complex conjugate ones appears in 
several regions of the spectrum. We have not been able to study whether this phenomenon persists 
as the system size L is increased, but it is a common feature in non-unitary models. 
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y 

FIG. 4. Same as Fig. [3] but for L = 4. 

In Figure 2] we present our data for the 7 dependence of the low energy spectrum of the L = 4 
mixed superspin chain. As 7 — > the spectrum (ordered by energy) consists of the E = 
singlet ground state, two degenerate octets, an 8-dimensional indecomposable and another pair of 
degenerate octets [4|. The observed splittings and appearance of levels with complex conjugate 
eigenvalues fit into the scheme discussed for L = 3 above. 



B. Analysis of the finite size spectrum antiferromagnetic regime 



Due to the identification (|3.8[) of the spectrum of the XXZ spin-1 chain within that of the 
superspin chain we already know part of the finite size scaling amplitudes - relative to the E = 
eigenstate - in the latter: 

A^-) = ^(2A ( ^( 7 )-i). (4.4) 

In translating the energies from the XXZ model we have used that as a consequence of (|3.7p the 
Fermi velocity remains the same, Vp 1 = vf = 27r/7, and therefore the scaling dimensions of the 
superspin chain are twice of those given in (|3.12p for the composite fields in the XXZ spin-1 chain. 
This is in agreement with the observations in Ref . |4 . 
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For 7 < 7r/2 the ground state is the unique state with Eq(L) = 0, hence the central charge of 
the model is 

J ^0 = 0, 0< 7 <^. (4.5) 



ttvf 2 

The lowest excitation corresponds to the conformal operator with scaling dimension xQ'®) = X c 
in the ir- twisted XXZ spin chain of odd length L (|3.15p . In the mixed chain this translates into a 
scaling dimension with scaling dimension 2X C — 1/4 = (ir — 2^)/4ir according to (|4.4p . 

In the superspin chain this state is realized by (L — l)/2 narrow 2-strings. As mentioned 
above, a narrow string may be viewed as degeneration of two strange 2-strings of opposite type. 
Excitations can be created by lifting this degeneracy into configurations with different number of 
the two possible types of strange 2-strings, i.e. type '+'-strings with A 1 - 1 ** = (A^ 2 ))* = A + 
and '-'-strings with AW = (A^ 2 ))* = A - 17/4. The narrow string state is described by the same 
number of ± strange strings, AiV = N+ — iV_ = 0. Similar as in Ref. U one can show that states 
with different but finite AiV have the same energy in the thermodynamic limit (note, however, that 
configurations with AiV even (odd) are only possible for L odd (even)). Computing the energies of 
these states we find a logarithmic fine structure on top of the level corresponding to the operator 
(n, m) = (1, 0), i.e. 

_J±_ AE (mix) = 1_z2lL + Kh , L)(AiV) 2 , (4.6) 

2lTVF 47T 

see Fig. [5] for 7 = 2-/r/7. We have also displayed the L-dependence of the AiV = level correspond- 
ing to the twisted spin chain excitation with (n,m) = (0, 1) for even L (|3.14p : this state leads to 
a scaling dimension j(tt + j)/(tt — 7) in the spectrum of the mixed superspin chain and is part of 
the lowest sZ(2| 1) indecomposable for even L as 7 — > 0. 

Starting from a non-linear sigma model with a supersphere as target space Ikhlef et at 
have proposed that this logarithmic fine structure in the finite size spectrum is the signature of a 
non-compact boson in the continuum theory. At intermediate scales, given by the size L of the 
lattice regularization, the renormalization of the coupling constant leads to an effective radius of 
the non-compact boson given as 

(Rnc? oc [ln(L/L )] 2 . (4.7) 

This, in turn, generates the logarithmic structure (|4.6|) observed in the numerical data: 
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FIG. 5. Evolution of the fine structure of the low energy spectrum of the mixed spin chain as a function 
of L for 7 = 2ir/7: circles denote energies scaling to the (n,m) = (1, —1) level of the XXZ spin chain with 
odd length, triangles denote the (0, 0) level of the XXZ chain. The dotted lines connecting to L = oo are 
rational function extrapolations of the numerical data. 

In Fig. O we present numerical estimates of the amplitude ^(7) for the lowest excitation based 
on the finite size spectrum for system sizes up to L = 4095. We find that the amplitude is well 
described by 

5 7T — 7 



A( 7 ) 



(4.9) 



2 7T + 7 

For 7 •> 1 the extrapolation of the numerical data is within 1% of this conjecture. For more 
convincing evidence, in particular at small values of the deformation parameter 7, one would need 
additional information about higher order corrections to (|4.8p and have to study system sizes which 
are out of reach for this method based on the numerical solution of the Bethe equations (|2.17p . 

For 7r/2 < 7 < 7r the set of levels ()4.6p has energies below the E = eigenstate of the mixed 
chain for L sufficiently large, see Fig. [71 As a consequence the model is in a different universality 
class with central charge determined by the finite size scaling behaviour of the lowest level in this 
set, i.e. 



c e ff(7) 



6L 

TTVp 



E {mix) = o 2 7 zJL 1 
7T ' 2 



< 7 < 7T . 



(4.10) 
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FIG. 6. Amplitude of the logarithmic fine structure extracted from comparison of the energies of corre- 
sponding states at sizes L = 2L' + 1 (2L' for L even) and L'. The line is the conjectured 7-dependence 
Eq. gU). 



Note that this state can be identified with the lowest level of the XXZ spin-1 chain for odd L only. 
For even L its realization in the mixed superspin is in terms of a strange string configuration with 
A./V = ±1. Immediately above this ground state there is a continuum of states with energy gaps 
vanishing as oc (AN) 2 / [L ln(L/Lo) 2 ] due to the presence of the non-compact boson: our finite 
size analysis indicates that these gaps show the same 7 dependence as in (|4.9|) given above. The 
first excitation from the XXZ spin-1 subset of the spectrum above this continuum is the E = 
state corresponding to a scaling dimension 



X 



27 — IT 
47T 



(4.11) 



V. FERROMAGNETIC REGIME OF THE MIXED CHAIN 



Again we use the spectral relation (|2.19p to discuss the properties of the mixed superspin chain 
in this regime in the interval < 7 < ir by changing the sign of the energy eigenvalues (|2.18p . 
As for the antiferromagnetic regime above we begin our analysis based on numerical data for the 
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FIG. 7. As Fig. [S]but for 7 = 2ir/3. The dotted line denotes the energy level at E = of the singlet state 
without any finite size scaling in the spectrum of the mixed chain for L both even and odd. 



spectrum of the L = 3 chain obtained by exact diagonalization of the superspin Hamiltonian, see 
Table [III for 7 = 2tt/7 < ir/3: according to our discussion of the disordered ferromagnetic regime 
of the XXZ spin-1 chain above the lowest state from the XXZ spin-1 part of the spectrum for this 
value of 7 is the completely polarized reference state with S3 = 3 and energy E = — 12cot(7/2) 
according to (|3.8p . In Table ITT1 we have also listed the level corresponding to the primary field 
(n,m) = (0, \) in the twisted XXZ model (pTTo]) . 



The numerical data reveal that below the reference state of the XXZ spin chain there are many 
states with lower energies in the spectrum of the mixed superspin chain in the charge-sectors 
b = 0, ±1/2 and with magnetization < S3 < L: the corresponding Bethe configurations consist 
of (N\, N2) roots with Ni — N2 = 0, ±1 and Im(X^) = n/2, i.e. (1, — )-strings. The ground state 
in the sector (L,L — 1) (together with the equivalent sector (L — 1,L)) is four-fold degenerate. 
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TABLE II. Low energy states of the L 
the identified Bethe configurations. 



3 superspin chain in the ferromagnetic regime for 7 = 27r/7 and 



{Ni,Ni) 


Energy E 


degeneracy 


Bethe roots 


(3,2) 


-27.8683 


4 


AW = 


{±1.1935 H 


- wr/2, itt/2}, A< 2 ) = {±0.4797 + iir/2} 


(3,3) 


-27.7508 


2 


aw = 


{-0.7367 H 


- i7r/2, 0.2722 + wr/2, 00} = -A< 2 ) 


(2,2) 


-27.4205 


4 


aw = 


{ — 1.0384 H 


- wr/2, 0.2015 + in/2} = -A( 2 ) 


(2,1) 


-27.1935 


4 


AW = 


{±0.5829 H 


-wr/2}, A< 2 ) = {ztt/2} 


(1,1) 


-26.4469 


4 


AW = 


{—0.4447 H 


-i7r/2} = -A< 2 > 


(1,0) 


-25.8814 


4 


AW = 


W2} 




(0,0) 


-24.9183 


2 


(XXZ 


pseudo vacuum) 


(3,2) 


-24.8425 


8 


A(D = 


{ — 1.0326 H 


- ztt/2, 0.3894 + wr/2, 0.6432} 








AW = 


{—0.2736 H 


-wr/2, 0.7440} 


(3,3) 


-24.7595 


4 


AW = 


{—0.5477 H 


-wr/2, 1.0369, 00} 








aw = 


{-00, -0.0189 + i7r/2, 0.5909} 


(2,2) 


-24.6243 


8 


AW = 


{-0.8222, 


-0.2226 + i7r/2} 








AW = 


{-0.7559, 0.8563 + i7r/2} 


(3,3) 


-23.6547 


2 


AW = 


{0.6023 ± i0.6800, 0.7226} = A< 2 > (XXZ) 



A. Thermodynamic limit 



Based on this observation we shall now study this ground state in the thermodynamic limit 



(5.1) 



L — > 00. To make further progress it is convenient to rewrite the rapidities as, 

where fij € R for a = 1,2. Now, by substituting Eq. (|5.1|) in the Bethe ansatz equations (|2.17|) 

(k) 

and by taking their logarithms we find that the resulting relations for ^l - are, 



N 2 



L$(^, 7 - vr) = 2„Q?> + £ - , 7) , j = l,..., N, 

k=l 
Ni 

L*{nf ,7 - vr) = 27rQf + J] $( M f - $\ 7) , i = 1, • • • , N 2 

fc=l 

»(«) 



(5.2) 



where $(x, 7 ) = 2 arctan (tanh(x) cot(7/2)). The numbers define the many possible branches 
of the logarithm. They have to be chosen integer or half-odd integer depending on the parities of 
N a according to the rule 



Qf = ^ mod 1, Qf = — mod ! . 



1(2) _ 

2 



(5.3) 
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For example the ground state in the sector N\ = L and N2 = L — 1 is described by the symmetric 
sequences 

) (D_ L + 1 „■ 



3 

(2) 
3 



J 



1,...,L, 



L 
2 



(5.4) 



3i 



J 



1,...,L-1. 



At this point we have the basic ingredients to compute some of the thermodynamic limit prop- 
erties. When L — > 00 the number of roots p,j tend towards a continuous distribution on the real 
axis whose density which we shall denote by p^ a '{p). The Bethe equations (|5.2p turn into coupled 
linear integral relations for the densities p^ (p) which can be solved by Fourier transform. The fact 
that Eqs. (|5.2[) are symmetric under the exchange of rapidities p^ o p? in the thermodynamic 
limit implies that p^\p) = p( 2 \p). The final result for such density is, 



/> (a V) 



I COS 


7T7 


cosh 


7T/i 






7T-7/2 


(tt - 7/2) cosh 


2irfi 


+ COS 


7T7 


71--7/2 


2(7r-7/2)J 



for a = 1,2. 



(5.5) 



Now from the expressions for the density p^ (//) and Eq. (|2.18p we can compute the ground state 
energy density eoo = Eq/L. By writing the infinite volume limit of Eq. (|2.18p in terms of its Fourier 
transform we find 



-4 cot (7/2) - 4 / do; 



sinh [077/2] cosh[cj7/4] 



for < 7 < 7T . 



(5.6) 



sinh[w/r/2] cosh[oj(27r — 7)/4] 

In addition, we have verified that the low-lying excited states about the ground state are gapless. 
As usual, these states can be obtained by inserting holes in the density distribution of fj,j by 
making alternative choices for Qj . This procedure is nowadays familiar to many integrable models 



solved by Bethe ansatz and for technical details see for example ([36438]). We find that the low- 
momenta dispersion relation among the energy e*- a ^ (/z) and momenta (p) , both measured from 



the ground state, has a relativistic behaviour 



-(mix) ( a )/ \ 



(5.7) 



The common slope at p^ a \p) = is the corresponding Fermi velocity of the excitations. It is 
determined by 



-(mix) 



(M) 



2vrp( a )(ii) 



2tt 



11=00 



2tt — 7 



(5.8) 



As expected from (|3.7p it coincides with the velocity vf of low energy excitations of the XXZ spin-1 
chain in the disordered ferromagnetic regime, see Section IIIIB1 
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B. Analysis of the finite size spectrum ferromagnetic regime 

From our investigation of the behaviour of the Bethe ansatz roots associated to the low-lying 
excitations we found that they can be well described in terms of real rapidities. The roots with fixed 
imaginary part at can easily be mapped on real roots by means the straightfoward shift (|5,1|) . 
We remark however that some excitations have the peculiar feature that some of the their roots 
have the real part located at infinity. This scenario suggests us that a first insight on the structure 
of the finite-size corrections can be obtained by applying the standard density root method, see for 
instance 3814421] . This technique explores the Bethe ansatz solution and it allows to compute the 
O (L~ 2 ) corrections to the densities of roots p^ k \n). This approach predicts that the finite-size 
corrections to the low-lying energies eigenvalues have the following form, 

+ o(L" 1 ), (5.9) 



2-rvvp 

E(L,j) - Leoo 



.1 + X mi ' m2 h) 



(5.10) 



L 

where the scaling dimensions X™^^ 2 (7) depend on the anisotropy 7 as 

In (|5.10p the integers n\ and ni are related to the number of roots at each level of the Bethe 
equations by N\ = L — n\ and N2 = L — Therefore (rti ± n-2)/2 are associated to the conserved 
U(l) spin S3 and charge b of the model: the scaling dimensions ()5.10j) exhibit exact spin charge 
separation in the low energy effective theory. The corresponding excitations of the model are free 
bosons with compactification radii i? 2 ~ 7/27r and ~ (1 — 7/271"), usually denominated spinon 
and holon modest The indices mi (m-2) are related to macroscopic momentum of the excitation 
due to backscattering processes on the first (second) level of the Bethe ansatz and they are usually 
called vortex excitations. As a consequence of (|5.3p they take integer (half-odd integer) values 
depending on the parity of N\ + N2 leading to the following constraint connecting spinon and 
vortex numbers: 

• for 77-1 ± 77-2 odd —7- 777,1, 7772 = 0, ±1, ±2, .. . 

! 3 5 (5.H) 

• for 77i ± n 2 even ->■ mi, 777 2 = ±-, ±-, ±-, . . . . 

Let us now investigate the validity of the formulae (|5.10p for the conformal dimensions together 

with the selection rule (|5.1ip . In order to do that we have solved the Bethe ansatz equations for a 

1 A similar observation has been made in one of the critical phases of a Temperley Lieb model with staggered spectral 
parameters where the effective field theory consists of a compact boson and two Majorana fermions [4^. In this 
model, however, the low energy degrees of freedom cannot be related to (7(1) charges of the microscopic model. 
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TABLE III. Finite size sequences 15.121 of the anomalous dimension X^(j) for 7 = 7r/6, 7t/3, 7r/2, 2-7r/3 
from the Bethe ansatz. The expected exact conformal dimension is X®'®^) = 4. 



^°n(7) 

1,U V / / 


7T 

6 


7T 

3 


7T 

2 


2tt 
3 


A 

4: 


f) 9fl4874fi7 




94&1 087^ 


D 94^41 37 


Q 
O 


PI 9' : t' : !7 c in71 


94^91 ^04 

U . Zr 4: U £1 J- O U ^± 


O 94Q7&&97 




1 9 


U.^4;0-LC'4;4;0 


9484^080 


94QQ991 


D 9407^8^ 


16 


24643218 


24913601 


24995927 


24987579 


20 


0.24780188 


0.24945223 


0.24997482 


0.24992223 


24 


0.24850243 


0.24962151 


0.24998285 


0.24994660 


28 


0.24891117 


0.24972275 


0.24998754 


0.24996104 


32 


0.24917161 


0.24978814 


0.24999053 


0.24997031 


Extrap. 


0.2503(2) 


0.250003(1) 


0.250002(2) 


0.250001(2) 


Exact 


0.25 


0.25 


0.25 


0.25 



number of low- lying states up to L = 32. From the numerical data we compute the sequence 

X{L) = (E(L, 7 ) - Leoo) + i (5.12) 



which in the thermodynamic limit is expected to extrapolate to the dimensions (|5.1Up . 

In Table ITlTI we show the finite-size sequences (|5.12p for the ground state Eq(L, 7) in the case 
of various values of 7. The extrapolated value of the corresponding conformal dimension turns out 
to be independent of the anisotropy 7 whose value is in good accordance with the one predicted 
by Eq. ([5T0D for X°'°( 7 ) = X°; 1 °( 7 ) = \. From this result we find that the ground state energy 
scales as 



c 



irvp 
~6L 



+ o(L-i), 



(5.13) 



which leads us to conclude that the continuum limit of the superspin chain in the disordered 
ferromagnetic regime should be described by a conformally invariant theory with central charge 
c = —1. The respective anomalous dimensions X^^ 2 (t) °f the theory have to be measured from 
the ground state (|5. 13[) and therefore they should be given by, 



X™ 2 (7) 



X mi,m 2 ( \ _ t 



(5.14) 



We have also analyzed the corrections to scaling due to the presence of irrelevant operators in the 
lattice Hamiltonian 44]. For all values of 7 we found the leading corrections to the finite size 
estimate (j5. 12|) for Xq '^(7) to be of order L~ 2 arising from the conformal block of the identity 
operator which explains the good convergence of the extrapolation. 
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TABLE IV. Finite size sequences 15.121 of the anomalous dimension Xq' 2 (7) for 7 = ir/6, 7r/3, 7r/2, 2tt/3 

1 _ 1 

from the Bethe ansatz. The expected exact conformal dimension is Xq' q 2 (7) = it /(An — 27). 



1 1 

A 0,0 17J 


TT 

6 


TT 

3 


TT 

2 


Itt 
3 


4 


0.23109569 


0.28352514 


0.33170562 


0.36094436 


8 


0.25745093 


0.29695955 


0.33325340 


0.37254539 


12 


0.26633674 


0. 29875319 


0.33331810 


0.37397623 


lb 


0.zb938340 


U./99315 ID 


U.ooooz857 


0.37443588 


20 


0.27066931 


0.29956684 


0.33333139 


0.37464230 


24 


0.27132594 


0.29970094 


0.33333245 


0.37475283 


28 


0.27170870 


0.29978104 


0.33333283 


0.37481895 


32 


0.27211767 


0.29983274 


0.33333303 


0.37486165 


Extrap. 


0.27273(2) 


0.300004(2) 


0.3333332(2) 


0.375001(2) 


Exact 


0.272727 . . . 


0.3 


0.333333 • • ■ 


0.375 



We now turn our attention to the excited states in order to bring extra support to the proposal 
(|5.10p . (|5.1ip . The first excitation above the ground state (|5.13p occurs in the sector n% = n<i = 0. 
Our analysis of the L = 3 system (see Table [TTJ) indicates that among the corresponding Bethe 
roots there is a pair of rapidities (A 1 - 1 ), X^) which takes values (±oo,=Foo). The presence of such 
infinities leads to an effective scattering phase shift of exp[±i7] for the remaining roots. Due to 
this peculiarity we shall present explicitly the form of the Bethe equations for the finite roots 

L-2 



L*( M f , 7 - tt) = 2vQf +7 + E ~ $ > 7) > 



vr) = 2nQf - 7 + £ 



k=l 
L-l 



(2) 



1,...,L-1 



,L- 1. 



(5.15) 



fe=i 



For this state the numbers Qj are given by 



Qf 



L + l _ 
2 

L + l 



J- 



+ 3, 



3 = 1, • • 
j = l, 



,L- 1. 



L-l 



(5.16) 



In table pvp we present the finite-size estimates associated to the state described by Eqs. (|5.15p 

with (|5.16p . We observe that the the extrapolated values agree with the proposal (|5.10p . (|5.1ip 

_ 1 _i 

which predicts that the lowest conformal dimension is Xq' 2 (7) = 4 ( 2 J_ 7 ) • 

In Tables IVllVIIII we present the finite-size sequences for several low-lying excitations in other 
sectors. We observe that the extrapolated values corroborate the result predicted by (|5.1U|) . 
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TABLE V. Finite size sequences 15.121 of the anomalous dimension Xf[ 2 (7) for 7 = tt/6, 7t/3, it/2, 2tt/3 

1 _ 1 

from the Bethe ansatz. The expected exact conformal dimension is 2 (7) — ^ + 2(2^-7) • 



1 1 

x i,i (7) 


7T 

6 


7T 

3 


7T 

2 


27T 

3 


4 


0.31054918 


0.44669719 


0.58683657 


0.74133601 


8 


0.33974964 


0.46283122 


0.58433551 


0.71638581 


12 


0.34923933 


0.46507338 


0.58378663 


0.71189057 


16 


0.35248317 


0.46578873 


0.58358979 


0.71032924 


20 


0.35385561 


0.46610987 


0.58349789 


0.70960905 


24 


0.35455798 


0.46628188 


0.58344777 


0.70921857 


28 


0.35496794 


0.46638478 


0.58341748 


0.70898230 


32 


0.35522923 


0.46645125 


0.58339779 


0.70883087 


Extrap. 


0.35603(2) 


0.466665(2) 


0.5833334(1) 


0.708332(2) 


Exact 


0.356060 . . . 


0.466666 • • • 


0.583333 • • • 


0.708333 • • • 



TABLE VI. Finite size sequences 15.121 of the anomalous dimension Xl%(j) for 7 = tt/6, tt/3, tt/2, 2vr/3 
from the Bethe ansatz. The expected exact conformal dimension is ^2(7) = \ + 7/^ 



< 2 °(7) 


7T 

6 


7T 

3 


7T 

2 


2tt 
3 


4 


0.36660004 


0.56006032 


0.75496874 


0.95242531 


8 


0.39901997 


0.57907838 


0.75200345 


0.92788809 


12 


0.40929983 


0.58158419 


0.75094335 


0.92183973 


16 


0.41280124 


0.58237304 


0.75054074 


0.91961495 


20 


0.41428312 


0.58272535 


0.75034904 


0.91856536 


24 


0.41504200 


0.58291356 


0.75024350 


0.91798979 


28 


0.41548512 


0.58302599 


0.75017939 


0.91764084 


32 


0.41576759 


0.583098552 


0.75013758 


0.91741358 


Extrap. 


0.4166(2) 


0.583332(1) 


0.750001(2) 


0.916662(2) 


Exact 


0.41666 • • • 


0.58333 • • • 


0.75 


0.91666 ■■• 



Of rather special nature is the state considered in Table IVIIII this is the lowest excitation from 
the XXZ spin-1 part of the spectrum. As in the antiferromagnetic regime the finite size scaling of 
these states can be deducted using the spectral relation (|3,8p . Using the same reasoning as above 
the scaling dimensions X nm {^) of the spin-1 chain in the ferromagnetically disordered regime, 
Eq. (|3.16p . should appear doubled in the spectrum of scaling dimensions of the superspin chain. 
Comparison with (|5.1(jp shows that this is indeed true (remember that m in (|3.16|) takes half-odd 
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TABLE VII. Finite size sequences 15.121 of the anomalous dimension X\' Q 1 ( r y) for 7 = 7r/6, 7r/3, tt/2, 2-7r/3 
from the Bethe ansatz. The expected exact conformal dimension is Xj L 'J" 1 (7) = j + / 2 j- ) ■ 



^1,0 V 1/ 
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6 
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2tt 
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Z. / oyDoo4o 
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l.oyU/ IvZo 






i.yyyiyo7o 
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1.709y0(54 
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1 £0Q07777 
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1.0104 fuoy 


l.oooolUo 


1. / OZ14004 


1 fi 

1U 


i.oooyuu'iy 


1 • z ±OOZU t ±UZ 


1 ^84800^ 
1 .oo^ouyo 


1 73QQZLQ89 


20 


1.48454418 


1.47493397 


1.5843369 


1.74356556 


24 


1.44315270 


1.46748697 


1.5840527 


1.74553075 


28 


1.41720088 


1.46292463 


1.5838718 


1.74671607 


32 


1.39992651 


1.45993403 


1.5837506 


1.7474855 


Extrap. 


1.3404(3) 


1.45003(1) 


1.5834(2) 


1.75002(1) 


Exact 


1.340909 • • • 


1.45 


1.58333 • • • 


1.75 



1 1 

TABLE VIII. Finite size sequences 15.121 of the anomalous dimension X 2 2 (7) for 7 = 27r/7, 27r/5, 27t/3, 

1 i_ 

57r/6 from the Bethe ansatz. The expected exact conformal dimension is X^ 2 (7) = j-. 





2tt 
7 


2tt 
5 


2ir 
3 


5tt 
6 


4 


1.84589078 


1.29492470 


0.75018486 


0.54556892 


8 


1.76868129 


1.25924483 


0.75002911 


0.58757296 


12 


1.75782587 


1.25398066 


0.75001177 


0.59492733 


16 


1.75432209 


1.25248723 


0.75000640 


0.59722733 


20 


1.75274406 


1.25141093 


0.75000403 


0.59824616 


24 


1.75189751 


1.25097723 


0.75000277 


0.59878935 


28 


1.75139056 


1.25071683 


0.75000203 


0.59911369 


32 


1.75106290 


1.25054825 


0.75000154 


0.59932298 


Extrap. 


1.75002(2) 


1.25003(1) 


0.7500002(3) 


0.600003(1) 


Exact 


1.75 


1.25 


0.75 


0.6 



integer values for twist ip = tt): 

= 2X n , m = n 2 2. + m 2 ^ (5.17) 
is exactly the pure spinon part of the low energy spectrum of the mixed superspin chain. 

i i 

As an example, we have computed the finite size sequences (|5.12p for the dimension Xq'q (7) = 
^. This is the state discussed at the end of Section fill Bl its root configuration changes at 7 = ir/k, 
k = 2,3,... evolving into single narrow string of L roots on both levels of the Bethe ansatz for 
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7 < ir/L. 

In the isotropic limit, 7 — > 0, this state is an s/(2|l)-descendent of the (completely polarized) 
reference state of the superspin chain with charges (6, S3) = (0, L) and lies outside of the low energy 
part of the spectrum. Unlike the XXZ spin-1 chain the superspin chain remains conformal in the 
isotropic limit of the disordered ferromagnetic regime: the holon sector, i.e. states from ()5.9|) . (|5.10p 
with mi = —ni2, remains conformal for 7 = 0. Taking into account (|5.1ip their scaling dimensions 
are 



ni.no \ 1 



2 (2n + l) 2 + vn? if rt\ — 712 = 2n + 1 and m\ 



m, 



(5.18) 



77,2 + I (^m + I) 2 if n i ~ n 2 = 2n and m\ = m + \ 

with integer m. These are the conformal dimensions of the isotropic osp{2\2) spin chain 

We have also studied the finite size scaling behaviour of several states with configurations which, 
apart from (1, — ) strings, contain real roots and inifinite roots on one or both levels (the existence 
of such configurations in the low energy sector of the ferromagnetic regime is indicated by our small 
system analysis, see e.g. Table Hi]) . In all cases we have considered these levels corresponded to 
descendents of the primary conformal fields identified before, i.e. with scaling dimensions X™^]^ 2 +n 
with integer n. 



VI. CONCLUSION 



In this paper we have studied an integrable f7 g [sZ(2|l)] vertex model built from alternating 



Based on its 



fundamental and dual three-dimensional representations first introduced by Gade 
solution by means of the algebraic Bethe ansatz we have computed the properties of this model 
in the thermodynamic limit and analyzed the finite size scaling of the low energy spectrum. From 
the latter we conclude that the critical point with central charge c = of the undeformed model 
identified before |J| is stable against variation of the anisotropy 7: as long as 7 G [0, vr/2) the 
ground state energy vanishes exactly without any finite size effects. In the continuum limit we 
find that the model displays a continuous spectrum of exponents in the sector with J7 9 [sZ(2|l)]- 
charge 6 = 0. As in the isotropic sZ(2|l) superspin chain the lower edges (|4.4p of the continua can 
be identified with the scaling dimensions in an antiperiodically twisted spin-1 chain, in this case 
the Fateev-Zamolodchikov model. The continuous part of the conformal spectrum leads to a fine 
structure (|4.6p in the spectrum of the large but finite superspin chain. This fine structure can 
be explained as signature of the presence of a non-compact degree of freedom in the continuum 
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theory with coupling constant renormalized to an intermediate scale of the order of the length L 
of the superspin chain. The 7-dependence of the coupling constant has been identified based on 
our numerical solution of the Bethe equations. At 7 = tt/2 a level crossing occurs changing the 
ground state and leading to an effective central charge taking values < c e ff (7) < 3 depending on 
the deformation parameter 7 6 [tt/2, it). Again, this critical behaviour mirrors that of the XXZ 
spin-1 chain of odd length L in the antiferromagnetically disordered regime subject to antiperiodic 
twisted boundary conditions. 

The spectrum in the ferromagnetic regime it < 7 < 2tt (or, equivalently, that of the chain 
with opposite sign of the exchange constant for anisotropies 7 = 2tt — 7) of the superspin chain 
is completely different: our finite size scaling analysis indicates that it is the same as for the 
C/ 9 [osp(2|2)] spin chain with central charge c = — 1. It displays separation of spin and charge 
degrees of freedom in the low energy excitations with the spin part of the spectrum turning non- 
relativistic as 7 — > 2ir. Unlike in the antiferromagnetic regime there are no signs of a non-compact 
degree of freedom in the continuum limit: the zero charge sector of the low energy spectrum can 
be identified exactly with that of the Fateev-Zamolodchikov model in its disordered ferromagnetic 
phase. 

In summary we have presented a comprehensive study of the critical properties of the mixed 
C/ ? [sZ(2|l)] superspin chain. The appearance of non-compact degrees of freedom in the continuum 
limit of such lattice models has been shown to be stable against deformation although it is limited 
to the antiferromagnetic regime of the model. As for the staggered six-vertex model studied in 
Ref. |6| our computation of the corresponding coupling constant (|4.9p relies on the numerical solution 
of the Bethe equations and its derivation within an analytical approach remains an open problem. 
More general phases can be expected to be found when one considers mixed chains based on 
higher-dimensional representations of the superalgebra (and its deformation). We note, however, 
that already in the corresponding XXZ spin-S" chains this leads to a growing number of phases 



(unitary and even non-unitary) as the deformation parameter is varied 
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